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ABSTRACT _ 

This paper represents a hybrid energy system (HES) consisting of 
photovoltaic (PV), Solid Oxide Fuel Cell (SOFC), electrolyzer system and a 
storage tank. In this proposed system a fuel cell controller is used where a 
PID controller is utilized to control the flow of hydrogen (Eb) through the 
valve to the SOFC to fulfill the load requirement. In this model a supervisory 
controller is used to regulate the whole system according to load 
requirement. So, when PV power is more than load requirement, then PV 
power fulfill the required load demand as well as the extra power of PV is 
utilized to generate the hydrogen (Eb) by the help of electrolyzer, further this 
hydrogen is used as a fuel of SOFC. Also, in this proposed model different 
types of faults are considered and verified their effect on the load as well as 
in distribution network. The computer simulations are done for the purpose 
technology and verified its effectiveness. Hence the appropriate results 
are obtained. 
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1. INTRODUCTION 

Renewable energy sources (RESs) provides clean energy sources which is pollution free, also RESs 
are environmentally friendly. These renewable energy sources are cost free and continuous according to [1], 
[2]. H. Yang, et al. [3] describes about various RESs consists of wind energy sources (WESs), photovoltaic 
energy system (PV), wave energy sources and micro-turbines, these are used in hybrid energy systems 
(HESs) and also applicable in micro-grid system. These RES structures are not operated properly due to the 
various hazardous conditions (solar radiation, temperature, wind speed). Hence numerous solutions are found 
out to solve these problems. Hence storage systems are the best solution for the mitigation of the effect of 
solar, wind and other renewable sources. Storage systems are used to store the energy, also improves power 
quality. It is convenient to keep more than one storage system to minimize the variation of power as well as 
the storage systems provide high energy density [4], [5]. Hybrid Energy Systems (HESs) uses Battery Bank 
(BBs) and storage systems efficiently, but the life time durations of these batteries decrease according to its 
charging and discharging periods [6], [7]. Hence an alternative secondary energy sources is required which 
can replace the battery bank as well as increase the efficiency of Hybrid Energy Systems (HESs). So fuel 
cells are used now a days, these fuel cells are combine with electrolyzer and gives power supply continuously 
to the load given in [8], [9].Hence the hybrid energy system (HESs) consists of PV, wind energy (WT) and 
fuel cell (FC) combine with the electrolyzer gives the better energy management strategies and more 
effective in practice for the requirement of higher energy and the uses of DC loads are more promising in 
today’s world such as they are used in commercial residential and industrial buildings. So, these loads are 
DC dominated in the mean time [10]. Hence the AC grids are used to supply these DC loads. So, DC/DC and 
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AC/DC converters are used more in the nearer days. So, the cost of the system is increased and gives the 
additional losses. So, the system efficiency is reduced. A HES containing an individual AC grid as well as a 
DC grid, which can overcome the problem that the DC loads are connected to the DC grid given in [11]-[13] 
and AC are connected to the AC grid. Hence bidirectional converters are used in DC as well as 
AC grids [14]. 

The literature consists of a PV system, SOFC, electrolyzer and a storage tank which are simulated 
with each other, latter this combination connected the grid. Advantage of This technology gives better 
performances to the DC as well as AC grids. Hence hybrid energy system increases the quality of power 
management of the grid. This article concluded overall description as well as provides the information about 
the proposed technology. Detailed modeling and controlling strategy of the system equipments are given. 
Further the proposed technology is verified for its performances and effectiveness under the simulation 
consisting of several situations. 


2. EXPLANATION OF GRID CONNECTED HYBRID GENERATION SYSTEM 

Figure 1 shows the proposed model which consists of PV system with the maximum power point 
tracker (MPPT) interfacing with the DC/DC converter (boost converter) to increase the DC voltage. DC/AC 
converters (Inverters) are used to convert DC to AC, when DC voltage of PV with MPPT having boost 
converter (Vpv_boost) is applied to this. Energy coming from the sun depends according to the various 
weather conditions, so, the solar radiation is not constant in every period, it may be weak or strong or absent. 

So, the integration of SOFC in the HESs makes the HESs more sustainable. When the Power 
generation is greater than the required amount, then the electrolyzer come to the front to take excess 
generated power and converts the water in to hydrogen and oxygen by electrolysis process. Storage tank 
stores this hydrogen in compressed form and used as a fuel of SOFC. 

Hybrid Energy System (HES) develops a power management strategy, which can decrease the effect 
of daily or seasonal variations due to the different climatic conditions. Hence HES provides a better output 
power, which can fulfill the requirement of power. In this proposed technology different types of faults are 
given, hence the stability and performances of the HES is analyzed. 



Figure 1. Hybrid system connected to grid (proposed model) 


3. PHOTOVOLTAIC (PV) SYSTEM 

The PV system mathematical model is represented in [15] and the maximum amount of power can 
be withdrawn from PV array with the help of MPPT; then DC/DC converter (boost converter) are used to 
boost up the voltage level and connected to the DC/AC inverter. 

Then the PV system is designed by considering (1) with the help up MATLAB/SIMULINK 
software and Table 1 shows the specifications of the photovoltaic (PV) array. 
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Parameters are considered for designing of the PV system such as, Io is the reverse saturation current 
of PV model in [A], Isc is the short-circuit PV cell current in [A], Ipy/I P his the output current of PV model in 
[A], k is the boltzmann’s constant in [J/°K], a is the completion or ideality factor, q is the electron charge in 
[C], Rp in the parallel resistance of the PV model in [Q], Rs in the series resistance of PV model in [Q],Ns is 
the no. of cells are connected in series in a string of the PV model, Np is the no. of parallel strings, T is the 
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temperature in [K], Vpy is the terminal voltage of the PV model in volt in [V], Vmp is the maximum power of 
the PV cell in volt in [V], VOC is the open-circuit voltage of the PV cell in volt in [V] 


Table 1. Specifications of Photovoltaic (PV) ARRAY 


Specification 

Value 

No. of series connected PV cell 

96 

No. of Modules in parallel 

66 

No. of Modules in series 

05 

Open circuit voltage (Voc) 

64.2 V 

Short circuit current (I S c) 

5.96 A 

Maximum Power of PV cell 

305W 

Temperature(T) 

25 

Irradiance 

1000 


4. MODELLING AND SIMULATION OF ELECTROLYZER 

Water is decomposed into two parts, one is hydrogen (H 2 ) and another one is oxygen, by electrolysis 
process when electric current is flowing through the electrolyzer, which contains two individual electrodes 
[16]. The equation of electrolysis of water is represented below: 

H 2 0 (/) + electrical energy = H 2 (g) + ~0 2 (g) (2) 

Parameters are used to design electrolyzer model such as, i e is the current of an electrolyzer in [A], 
n c is the no. of electrolyzer cells connected in series, nPE is the hydrogen generation in moles per second in 
the [mol s -1 ], F is the faraday constant in the [C kmol -1 ], r|Fis the Faraday efficiency. 

Generation of hydrogen of an electrolyzer is directly proportional to the electric current, which is 
related to the Faraday’s law, is given below [16]. The electrolyzer model is simulated with the help of 
simulink by using (3) and (4). 


nH 2 


n F n c i e 
2 F 


( 3 ) 


The ratio between hydrogen generations by the help electrolyzer to the actual value is given as 
Faraday efficiency. 
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5. MODELLING AND SIMULATION OF STORAGE TANK 

Electrolyzer is used to generate the hydrogen by the help of electrolysis process; then this hydrogen 
is used as the input fuel of SOFC. So SOFC consumes the required amount hydrogen according to the power 
requirement. The rest amount of hydrogen is stored in the storage tank after the use of hydrogen as a fuel to 
SOFC, which is the difference of the generated hydrogen and the required hydrogen used by SOFC. 

The hydrogen storage tank stores the hydrogen in compressed manner either it is stored in liquid or 
gasses manner. From (5) is used for simulation and design of a hydrogen storage tank according to [16]. 


?b ~ Pbi 


NH 2 RT b 
Z MH 2 V b 


( 5 ) 


Parameters are used to design the Storage Tank such as, Z is the pressure compressibility factor, 
NH 2 is the amount of hydrogen stored in the storage tank in moles per second in the [kmol s -1 ], MH 2 is the 
Hydrogen molar mass in [kg kmol -1 ], Tb is operating temperature in [°K], Pbi is the storage tank pressure at 
initial stage in [Pa], Pb is the storage tank pressure in [Pa], R is the Rydberg/Universal gas constant in [J 
(kmol °K) -1 ], Vb is the storage tank volume in [m 3 ]. 


6. DETAILED MODELLING OF SOFC 

Parameters are used to develop a mathematical model of the solid oxide fuel cell (SOFC), which is 
shown below according to [17] such as, B and C is the activation voltage constants in SOFC system in [V], 

CV is the conversion factor in r Hydrogen ( kmoh ) 1 g j s the normal no load voltage in [V], E is the nernst 

L Methen (kmoh) ]’ & L J ’ 
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immediate voltage in [V], F is the Faraday’s constant in [C (kmol' 1 )], Ifc is the feedback current of SOFC in 
[A], ki is the gain of PI (proportional-integral), K an is the anode valve constant in [^/kmol kg(atm s) _1 ], K r is 
the modeling constant in [kmol SA' 1 ], Kh 2 o is the molar constant of water valve in [kmol (atm s' 1 )], Kh 2 is the 
molar constant of hydrogen valve in [kmol (atm s" 1 )], K 02 is the molar constant of oxygen valve in [kmol 
(atm s' 1 )], Mh 2 is the hydrogen molar mass in [kg kmol' 1 ], Ph 2 o is the water fractional pressure in [atm], Ph 2 
is the hydrogen fractional pressure in [atm], P 02 is the oxygen fractional pressure in [atm], No is the no. of 
fuel cells are in series in a stack, qm is the molar flow of hydrogen in [kmols' 1 ], q me thane is the methane flow 
rate in [kmols 1 ], qo 2 in the input molar flow of oxygen, [kmols -1 ], qn 2 is the input flow of hydrogen in 
[kmols 1 ], q^ is the hydrogen output flow in [kmols 1 ], q^ is the hydrogen flow amount that reacts in 
[Kmols 1 ], q^ e 2 q is the hydrogen flow amount that meets load change in [kmols 1 ], R is the rydberg (Universal) 
gas constant in [(latm) (kmol K) _1 ], Rint is the internal resistance of SOFC in [Q], T is the absolute 
temperature in [K], U is the utilization rate,V an is the volume of the anode in [m 3 ], V ce ii is the SOFC dc output 
voltage in [V], iq x 2 is the reformer time constants in [s], x 3 is the PI controller time constant in [s], x H2 in the 
hydrogen time constant, [s], x H20 is the water time constant in [s], t 02 is the oxygen time constant in [s], Z act 
is the activation over voltage in [V], Z 0 hmic is the ohmic over voltage in [V]. 

The ratio between the molar flows of hydrogen through the valve to the hydrogen partial pressure in 
a channel can be given as [17]. 


QH2 k CLTL |t 

- 7^2 - H 2 


( 6 ) 


Molar flow of hydrogen is calculated by taking the derivative of the partial pressure, which shows 
the perfect gas equation as given in [17]. 


dt VH2 ~ (Qh 2 Rtn Rhi) 


(7) 


Relationship between hydrogen flow and SOFC current determines the rate of reacted hydrogen 
flow according to [17]. 


Rh2 ~ 


N 0 J FC 
2 F 


= 2 KM 


( 8 ) 


The partial pressure of hydrogen can be determined by considering the Laplace transform of (6)-(8) 

as [16]. 


Ph2 ~ 


1 /k i 


H 2 
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Where, 


T H2 ~ 


Van 

K H2 RT 


( 10 ) 


In this manner oxygen partial pressure and water partial pressure can be determined. The output 
voltage of the fuel cell can be obtained by considering the summation of three factors; these are activation 
overvoltage of the anode and cathode, Nernst potential and ohmic overvoltage of internal resistance 
as given in [17]. 


^cell — F + q act + q 0 hmic 

Where, r\ act = -Bln(CI PC ) 

And H ohmic — ^ Ifc 

So, Nernst instantaneous voltage can be determined as given in [17]. 
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The hydrogen acts as a fuel of SOFC power plant which can able to fulfill the power demand 
according to its requirement. Hence the supply of hydrogen is continuous for the stack operation of SOFC. 
The model of the reformer can be expressed as [17]. 


Qmethanol t 1 t 2 s 2 +(t 1 +t 2 )s+1 

A PID controller is used to control the flow of hydrogen according to the requirement of output 
power of SOFC system. SOFC output current is taken as a input of the SOFC to obtain the feedback control 
[17], given by 


req 

Qh2 


NqIfc 
2 FU 


(16) 


The available hydrogen of the reformer is used to control the flow of methane by the help of PI 
controller [17], which can be determined as shown below. 

^methane = (*1 + “ Qhz) (17) 

By considering the above parameters a 5 KW SOFC model is obtained as given in Table 2. A SOFC 
brief model is represented in simulink as given in [18]. Three 5 KW SOFC are used in this proposed system. 


Table 2. Specifications of SOFC to design the mathematical model 


Specifications 

Value 

Activation voltage constant (B) 

0.04777 [A* 1 ] 

Activation voltage constant (C) 

0.0136 [V] 

Conversion factor (CV) 

2 

Faraday’s constant (F) 

96 484 600 [Ckmol 1 ] 

FC system internal resistance (R int ) 

0.26664 [Q] 

FC absolute temperature (T) 

343 [K] 

Hydrogen time constant (x H2 ) 

3.37 [s] 

Hydrogen valve constant (Kh 2 ) 

4.22 xlO' 5 [kmol (s atm)' 1 ] 

Hydrogen-oxygen flow ratio(r H 2 o) 

1.168 

Kr constant =N 0 /4F 

2.2802 x 10' 7 [kmol (s A)' 1 

Line reactance (X) 

0:0809 [Q] 

Methane reference signal (Q me thref) 

0.000015 [kmol s' 1 ] 

No load voltage (E 0 ) 

1.18 [V] 

Number of cells (N 0 ) 

88 

Number of stacks (N s ) 

1 

Oxygen time constant (x 02 ) 

6.74 [s] 

Oxygen valve constant, (k 0 2 ) 

2.11 xlO' 5 [kmol (s atm)' 1 ] 

PI gain constant (ki) 

0.25 

Reformer time constant (x x x 2 ,x 3 ) 

15 [s] 

Universal gas constant (R) 

8314.47 [J kmolK' 1 ] 

Utilization factor (U) 

0.8 

Water time constant (x H2 o) 

18.418 [s] 

Water valve constant (Kh 2 o) 

7.716 10' 6 [kmol (s atm) 1 ] 


7. INTERNAL STRUCTURE OF HYBRID ENERGY SYSTEM 

Figure 2 gives the information about the block diagram of Hybrid Energy System (HES), which is 
represented in Figure 1. By using the MPPT in the PV system maximum amount of power can be tracked. 
Then the PV power is used to manage the power demand of the load as well as the rest amount of PV power 
is utilized to convert the water in to hydrogen and oxygen with the help of electrolyzer by the electrolysis 
process and this generated hydrogen is used as the input fuel of the SOFC shown in Figure 2(b). The DC/DC 
converter is connected to PV system with MPPT to increase the voltage of PV system (Vpv_boost). This 
DC/DC converter is triggered by the gate pulse which is generated from the MPPT as shown in Figure 2(a). 
Again, another DC/DC converter (boost converter) is connected to the SOFC to increase the voltage level of 
the SOFC voltage as shown in Figure 2(b). 

The PV system and SOFC boost converter are connected to common DC bus. Hence the DC voltage 
of the common DC bus is supplied to DC/AC converter (Inverter). So, the DC voltage is converted to AC 
voltage by the help of DC/AC converter shown in Figure 2(a). Then the output of the DC/AC converter is 
connected to grid. To synchronize the grid with hybrid energy system (HES) a voltage source inverter (VSC) 
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is used and DC/AC converter is triggered by the gate pulse provided by the VSC controller as shown 
in Figure 2(a). 

VSC controller is used to synchronize the HES voltage to the grid voltage. Hence grid voltage and 
frequency are utilized for the DC/AC converter operation. VSC is the combination of two cascaded control 
loop, one is voltage controller of DC link and another one is internal current loop. 

a. Voltage controller ofDC-link 

The desired value of active power is controlled by DC-link voltage. The id* is the output current of DC- 
link PI controller and it is used as the reference current or active current controller shown in Figure 2(c). 

b. Internal current loop 

The d->q control (Synchronous Reference Frame) is used for transformation of grid voltage (V a , Vb, 
V c ) and grid current (i a , ib, i c ) to dq reference frame which rotates synchronously along with grid voltage. 
With the help of q controller, the DC values are obtained by controlling the variables. So, for the 
designing and controlling purposes DC values are easier than AC values. Phase angle of the grid voltage is 
derived by using the PLL (Phase Lock Loop). PLL also used for transformation of abc^dq and also used for 
synchronization of grid voltage and controlled current. The phase angle is derived from the grid voltage by 
using the PLL technique. The value of reference reactive current (i q *) is taken as zero as shown in Figure 
2(c). So, power factor (pf) is maintained as one (Unity). The output voltages (Vd*, V q *) is provided by the PI 
controller and these output voltages are used as the input to the PWM generator, which produce pulses and 
given to the DC/AC converter to trigger [19]. 

From the Figure 2(b), a DC/DC converter (buck converter) is connected to the electrolyzer. So, the 
extra amount of power of the PV system is send to the electrolyzer through this buck converter. From the 
electrolyzer hydrogen and oxygen is extracted from the water and the hydrogen is used as the input fuel of 
SOFC. The hydrogen requirement of the SOFC depends upon the load requirement. Storage tank stores the 
rest amount of hydrogen, which is the difference between the generated hydrogen of an electrolyzer and 
required hydrogen of the SOFC. 

The block diagram of the fuel cell controller is shown in Figure 2(d). Voltage (Vsofc_boost) of the 
boost converter, which is connected to the SOFC, is applied to the PID controller. Then the valve is regulated 
by the PID controller according to the load requirement. So, the valve is decided to take the required amount 
of hydrogen through it and supplied to the SOFC. PID controller regulates the valve according to the load. 
The values are decided by the help of self tune mode of the PID controller, which are applied in PID 
controller (P=0.00013317, I=3.3925227e-05, D=-6.254254172870e-05). 



(a) Block diagram of hybrid energy system (HES) ( c ) Block diagram representation of 

voltage source inverter (VSC) 





(b) Block diagram representation of hydrogen flow from 
electrolyzer to SOFC 


(d) Block diagram representation of fuel 
cell controller 


Figure 2. A detailed structure of HES 


8. FAULT ANALYSIS OF THE PROPOSED SYSTEM 

Various types of fault such as symmetrical fault (Line to line to line fault or LLL fault, Line to line 
to line to ground fault or LLLG fault) and unsymmetrical fault (Line to ground fault or LG fault, Line to line 
fault or LL fault, Line to line to ground fault or LLG fault) occurs in the distribution network or at the load 
side which affects the supply voltage, current also affects the system performances according to [20], [21]. 
Hence the proposed system is verified under these various fault conditions. 
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So, fault analysis is performed on hybrid energy system connected to grid in MATLAB/Simulink 
environment as shown in Figure 1, hence the behavior of the proposed system is analyzed during the 
healthier condition and also the condition when fault occurs in the proposed system. In this proposed system 
both symmetrical as well as unsymmetrical fault are given and their effects are analyzed. 

9. OPERATIONAL CONTROL STRATEGY OF HYBRID ENERGY SYSTEM 

Operational control strategy (OCS) or supervisory control strategy (SCS) of the hybrid energy 
system is necessary to regulate the hybrid energy system as well as the load demand according to their 
effectiveness. Hence operational control strategy contains the four possible ways to regulate both the load 
demand and OCS. 

a. When Load demand (Pload) is smaller than the PV power (Ppv) [Ppv < Pload], at that time the required 
amount of PV power is used to fulfill the required load. Then the extra PV power is used for electrolysis 
process of the electrolyzer, where two basic elements such as hydrogen and oxygen are derived from the 
water. Hence the generated hydrogen is kept in a storage tank; later this hydrogen is used as a fuel of 
SOFC by flowing through the valve. 

b. When Load demand (Pload) is greater than PV power (Ppv) [Pload > Ppv], at that time to fulfill the 
required amount load, there is no sufficient PV power is available. So, the extra amount power is 
needed for the required load, which is provided by the SOFC. 

c. When the Load demand (Pload) is equal to the PV power (Ppv) [Pload = Ppv], at that time the load 
consumes the total available PV power. Hence the electrolyzer is not supplied from PV system. 

d. When generated PV power (Ppv) is zero [Ppv = 0], at that time, there is no supply available from the 
PV system. Hence SOFC is used to fulfill the load requirement. 

In this manner the four possible ways are find out for the operational control strategy to regulate the 
whole system. 


10. SIMULATION RESULTS AND DISCUSSIONS 

The proposed system (PV/SOFC/Electrolyzer/Storage tank) is designed by using the software 
MATLAB/Simulink, also different simulation is done and verified according to the system performances. 

Figure 3 shows the hydrogen generation and flow of hydrogen to the SOFC, also its flow is 
controlled by the PID controller through the valve. So, Figure 3(a) shows the hydrogen generation 
by the electrolyzer. Figure 3(b) shows the requirement of hydrogen flow to the SOFC according to 
the Load demand through the valve, which is controlled by the PID controller. Hence the PID 
controller output is shown in Figure 3(c). The output voltage of SOFC (Vsofc) is applied to DC/DC 
converter (Boost converter) as shown in Figure 2(b), where the voltage Vsofc is increased to 
Vsofc_boost (Voltage of boost converter connected to SOFC to increase the voltage Vsofc) as 
shown in Figure 3(d). 



Time (t) -> 

(a) H 2 generation by electrolyzer 


A 

O 

j 

a 


Time (t) 

(c) PID controller output 


Time (t) 

(b) H 2 flow to the SOFC. 


t 

8 

1 

5 

B 


Time (t) -> 

(d) Output voltage of boost converter connected 
to SOFC 


Figure 3. H 2 flow to SOFC and boost voltage of DC-DC converter connected to SOFC 
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Figure 4 shows the voltage and current wave form of the hybrid energy system (HES) on the 
healthier condition. But Figure 5 shows the voltage and current waveform of HES on L-G fault condition. 
Before fault occurs, the voltages of HES are remained constant but on the occurrence of L-G fault at the 
duration of 0.3 to 0.4 sec, then the voltage of faulted phase is equal to zero and other phases are in normal 
condition. When the fault is removed at 0.4 sec, the voltage of faulted phase is settled to its normal value 
after 0.4 sec. During the fault the currents of the faulted phases are increased. When the fault is removed; 
then the currents come to their original values. In this manner other faults of HES can be determined. 
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Figure 4. Voltage and current waveforms of hybrid 
energy system (HES) 


Figure 5. Voltage and current waveforms of HES 
during L-G fault 


Figure 6 gives the waveform of voltage, current and power of 30 load having the capacity of active 
power (P) =8000W, inductive power (Ql) =2000VAR and capacitive power (Qc) =100VAR. 

Figure 7 shows the symmetrical fault as well as unsymmetrical fault. From these Figures it can be 
seen before the fault occurs the voltage, current and power wave forms are with their appropriate values but 
when the fault occurs having the duration between 0.3 sec to 0.4 sec at that time the voltage, current and 
power wave forms of the faulted phase are equal to zero. After the fault is removed at the time of 0.4 sec, 
then the corresponding voltage, current and power wave forms. In this manner other faults of 30 load 
can be determined. 
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Figure 6. Voltage, current and power waveforms 
of 3-phase load 



Tune (t) 

Figure 7. Voltage, current and power waveforms 
of 3-0 load under L-G fault 


11. CONCLUSION 

From this paper a detailed model of hybrid energy system combination of PV/FC/Electrolyzer and 
storage tank is determined. The system configuration and the characteristics of the different components are 
verified as well as the power management strategy and overall control strategies of the system are given. In 
this proposed system the PV cell is the main generating source and used to fulfill the load requirement. 
Excess power of PV is used for the production of hydrogen from the electrolyzer. When the PV power is not 
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sufficient to fulfill the load requirement at that time SOFC is used to fulfill the required load by taking the 
hydrogen as a fuel; also, the amount of fuel flow through the valve is consumed by the SOFC is determined 
by the PID controller. 

Different types of fault analysis are considered on the proposed system and verified the effects due 
to these faults. Again, the performances of the system are verified by the help of simulations at different 
conditions. Further, the obtained results are conformed and provide the satisfactory performances. A better 
power management as well as control strategy are used for fulfillment of the the load requirement under the 
various critical and unfavorable situations. 

In this proposed system the hybrid energy system is connected to the grid with the help of VSC 
which gives the better results as well as better performances 
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